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A  Reynolds  stress  model  for  turbulent  boundary  layers  is  used  to 
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although  no  quantitative  explanation  is  offered  for  the  reduction  of  heat 
transfer  with  increasing  roughness  observed  by  Holden  at  M^  =  9.4. 

The  present  calculations  indicate  that  the  mean  velocity  is  approxi¬ 
mately  uniform  over  much  of  the  height  range  below  the  tops  of  the  elements, 
y  _<  k.  With  this  constant  "roughness  velocity,"  it  is  simple  to  estimate 
the  form  drag  on  the  elements.  This  roughness  velocity  has  been  investi¬ 
gated  by  systematically  exercising  the  present  model  over  ranges  of  poten¬ 
tial  parameters.  The  roughness  velocity  is  found  to  be  primarily  a  func¬ 
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I .  INTRODUCTION 


Surface  roughness  plays  an  important  role  in  turbulent  boundary 
layer'  skin  friction  and  heat  transfer  for  many  high-speed  flight  appli¬ 
cations.  Although  the  general  nature  of  roughness  effects  for  typical 
types  of  "sand  grain"  roughness  has  been  known  for  many  years,  dating 
back  to  the  classic  study  by  Nikuradse,^"  modern  composite  materials  in¬ 
troduce  a  different  character  of  roughness.  By  roughness  character,  we 
mean  the  shape,  spacing  and  perhaps  the  distribution  of  roughness  heights. 
With  a  woven  composite  material,  for  example,  the  exposed  fibers  would 
be  approximately  cylindrical  in  shape,  in  contrast  to  the  more  nearly  hemi¬ 
spherical  or  pyramidal  shape  of  conventional  roughness  elements. 

Previous  studies  of  roughness  character  have  been  somewhat  limited 

experimentally,  and  the  only  theoretical  investigations  have  been  quite 
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empirical.  Schlichting  measured  the  drag  due  to  various  element  shapes 
(spheres,  spherical  segments,  cones)  at  several  relative  spacings  on  the 
side  wall  of  a  water  channel.  There  are  several  reported  experiments 
involving  two-dimensional  roughness  elements  (machined  grooves  normal  to 
the  flow  direction).  Some  years  ago ,  Bettermann^  correlated  the  available 
data  to  obtain  k^/k,  the  ratio  of  the  effective  sand  grain  roughness  height 
to  the  actual  roughness  height,  as  a  function  of  roughness  shape  and  spac¬ 
ing.  Of  course,  the  effective  sand  grain  roughness  can  be  used  in  Niku- 
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radse's  results  to  predict  the  skin  friction  increase.  Dvorak  has  ap¬ 
plied  Betterman's  data  to  practical  heating  applications.  Unfortunately, 
much  of  the  data  that  had  previously  been  correlated  were  obtained  on 
two-dimensional  roughness  patterns.  One  might  expect  a  difference  in 
the  nature  of  the  flow  over  2-D  versus  3-D  roughness  .  For  example,  with 

2-D  roughness,  the  flow  would  be  more  likely  to  separate,  resulting  in 

.  ,  <> 

a  cavity  flow  in  the  grooves  between  the  elements.  The  three-dimensional 
case  is  of  far  greater  practical  importance,  and  the  2-D  type  of  roughness 
will  not  be  considered  further  here. 


A  second  important  issue  concerns  the  nature  of  roughness  effects 
in  supersonic  flow  conditions.  Density  variations  could  alter  the  extent 
of  roughness-induced  augmentation  of  friction  or  heating,  and  the  character 
of  the  flow  about  the  elements  might  change  if  supersonic  conditions  pre¬ 
vail  locally.  Only  recently  have  measurements  started  to  become  available 
on  supersonic  and  hypersonic  rough-wall  boundary  layers,  .offering  an  oppor¬ 
tunity  to  examine  the  Mach  number  issue  in  a  preliminary  fashion. 

In  this  study  we  employ  a  Reynolds  stress  turbulent  boundary  layer 
model  which  specifically  accounts  for  roughness  effects.  Roughness  is 
represented  by  distributed  sources  and  sinks  in  the  various  governing 
equations.  The  most  important  term  is  a  sink  term  in  the  mean  momentum 
equation  representing  form  drag  on  the  roughness  elements.  In  previous 
studies,'*  the  approach  was  developed  and  compared  against  subsonic  rough 
wall  boundary  layer  measurements.  The  present  objective  is  to  apply  the 
theory  to  variations  in  roughness  character,  as  well  as  to  supersonic 
conditions.  A  model  extension  to  treat  closely  packed  roughness  elements 
will  be  described.  Calculations  will  be  compared  with  a  variety  of  data 
to  establish  confidence  in  the  model;  then,  the  model  will  be  exercised 
systematically  over  a  range  of  parameters  to  develop  a  scaling  law  for 
roughness  character  that  is  simple  and  quite  different  from  existing  approaches. 


II.  ROUGH  WALL  TURBULENCE  MODEL 


The  turbulence  model  used  here  is  one  in  which  closure  approxima¬ 
tions  are  applied  at  second  order.  With  the  exception  of  the  treatment 
of  roughness,  the  formulation  is  somewhat  standard  at  this  time,  and  has 
been  successfully, applied  to  a  variety  of  smooth  wall  boundary  layer  and 
free  shear  flows.  The  model  accounts  for  both  mean  and, fluctuating  velo¬ 
cities  and  temperatures.  The  dependent  velocity  variables  are  the  mean 


velocity  vector  U^ ,  the  Reynolds  stress  tensor  u_  u^ ,  and  the  isotropic 
dissipation  rate  4>.  The  analogous  thermal  variables  (temperature  or, 
more  precisely,  enthalpy  h)  are  the  mean  enthalpy  h,  the  mean  square  fluc- 


.,2 


tuating  enthalpy  h'  ,  and  the  Reynolds  heat  flux  vector  u.'h' .  Under 


the  boundary  layer  approximation,  this  set  of  variables  reduces  to  U, 

2  2  2  2 

V,  u'  »  v'  ,  w'  ,  u'v',  4> ,  h,  h’  ,  u'h',  and  v'h'.  The  development  of 
the  governing  equations  and  the  required  closure  approximations  are  de¬ 
scribed  in  Ref.  5  and  will  not  be  repeated  here.  The  actual  equations 
are  given  in  the  Appendix. 


The  effect  of  roughness  is  described  by  distributed  source  or 
sink  terms  in  the  appropriate  equations.  As  already  noted,  only  distri¬ 
buted  roughness  is  considered  here,  and  we  make  the  fundamental  assumption 
that  the  flow  around  individual  elements  is  attached  to  the  elements. 

For  two-dimensional  roughness,  the  flow  might  be  treated  more  appropriately 
as  cavity  flows  between  the  elements.  In  the  present  model,  roughness 
elements  provide  a  distributed  sink  (due  to  drag)  for  mean  momentum,  and 
distributed  sources  for  mean  turbulent  kinetic  energy  and  dissipation. 

We  idealize  the  rough  surface  as  being  made  up  of  identical  elements  (al¬ 
though  the  extension  to  a  size  distribution  is  feasible) .  The  bottom  of 
the  elements  corresponds  to  y  =  0.  Let  k  be  the  element  height,  D (y)  be 
the  element  diameter  at  height  y  (for  0  _<  y  _<  k)  ,  and  Z  be  the  average" 
center- to-center  element  spacing.  The  functional  form  of  the  diameter 
D(y)  is  easily  prescribed  for  simple  shapes  such  as  cones  or  hemispheres. 

As  discussed  in  Ref.  5  and  in  the  Appendix,  the  drag  on  the  elements  per 
unit  volume  is  the  appropriate  sink  term  in  the  equation  for  the  mean 
velocity: 


(1) 


Ru  =  "  I  p  °2  cd  D(y)/i2  • 


A  drag  coefficient  value  of  CD  =  0.6  is  roughly  appropriate  for  elements 
such  as  cones  or  hemispheres.  The  source  terms  for  kinetic  energy  and 
dissipation,  which  are  less  important  numerically,  are  given  in  the  Ap¬ 
pendix.  Except  in  the  Stokes  flow  regime,  heat  transfer  to  an  element 
should  be  small.  Therefore,  the  only  roughness  term  appearing  in  the 
thermal  equations  is  a  source  for  the  mean  static  enthalpy.  This  term 
is  simply  constructed  so  that,  in  combination  with  Eq.  (1) ,  form  drag  does 
not  alter  total  enthalpy: 


Rjj  =  +  \  P  U3  CD  D  (y)  /£2 . 


In  the  approach  that  we  have  just  outlined,  roughness  elements 
are  assumed  to  occupy  no  physical  space.  This  assumption  becomes  pro¬ 
gressively  worse  as  the  roughness  density  increases.  Accordingly,  the 
model  has  been  extended  to  account  for  the  blockage  effect  of  the  rough¬ 
ness  elements.  This  is  done  in  a  simple  manner.  At  a  given  height  y, 
the  fraction  of  the  flow  area  in  the  x  direction,  that  is  open  to  the  flow, 
is  1  -  D(y)/£;  terms  that  act  in  the  streamwise  direction,  such  as  the 
convective  operator  pu  3/3x,  are  multiplied  by  this  factor.  Terms  that 

act  on  a  surface  area  whose  normal  is  in  the  y  direction,  or  that  act  on 

2  2 

a  unit  volume,  should  be  multiplied  by  l  -  ttd  /4£  .  However,  the  distri¬ 
buted  roughness  source  or  sink  terms  are  already  based  on  the  total  volume, 

rather  than  available  flow  volume,  and  need  no  such  factor.  If  the  entire 

2  2 

equation  is  divided  by  B(y)  =  1  -  ttd  /4£  ,  a  relatively  simple  change  occurs. 
For  example,  the  mean  momentum  equation  becomes 


, ,  .  3p  ,  1  3  „  3u 
f(y)  3x  B  3y  BU  3y 


h  (puV)  -  ip°2  CD ~^2 B_1 


where 


f  (y) 


1  -  D/l 
1  -  7TD2/4 i2 


(4) 


The  function  f (y)  contains  the  main  effect  of  blockage  and  may  be 
handled  by  merely  redefining  the  standard  stream  function  which  is  used 
to  eliminate  the  normal  velocity. 


||  -  f(y)P 


(5) 


Note  that  if  the  elements  are  packed  so  tightly  that  they  are  touching 
over  some  range  of  y,  then  D  =  £  and  f(y)  =  0.  The  stream  function  would 
be  forced  to  be  independent  of  y  over  that  range  (from  Eq.  (5) ) j  the  velo¬ 
city  would  remain  zero  from  the  bottom  of  the  elements  up  to  the  height 
where  D  <  £  and  flow  is  unblocked.  Of  course,  common  sense  would  dictate 
redefining  y  =  0  as  the  lowest  point  where  the  flow  is  unblocked.  How¬ 
ever,  the  model  does  yield  the  limiting  result  that  U  =  0  if  there  is  no 
space  between  the  elements. 

A  major  advantage  of  the  present  Reynolds  stress  model  is  that 
solutions  are  obtained  for  both  velocity  and  thermal  variables.  Heat  trans¬ 
fer  is  obtained  directly,  without  invoking  a  Reynolds  analogy.  Finite  dif¬ 
ference  solutions  are  obtained  using  the  obvious  boundary  conditions:  fluc¬ 
tuating  quantities  are  zero  at  the  solid  wall  and  in  the  free  stream.  It 

is  important  to  note  that  the  boundary  conditions  at  the  wall  are  not  treated 
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as  empirical  functions  of  roughness  as  has  been  done  in  other  approaches. 

For  numerical  solutions,  the  equations  are  first  transformed  to  the  stream 
function  coordinate,  guaranteeing  mass  conservation  and  eliminating  the 
normal  velocity  V.  The  transverse  coordinate  is  normalized  by  the  edge 
value  of  the  stream  function,  so  that  additional  mesh  points  need  not  be 
carried  in  the  free  stream  to  allow  for  boundary  layer  growth.  For  proper 
resolution  of  the  region  near  the  wall,  a  linear  mesh  in  the  logarithm 
of  the  stream  function  is  used.  The  finite-difference  equations  are  solved 
with  a  block  tridiagonal  Newton-Raphson  technique. 


III.  COMPARISON  WITH  EXPERIMENTS  ON  ROUGHNESS  CHARACTER 


The  first  measurements  to  be  examined  are  those  of  Schlichting, 
which  were  obtained  in  a  4  cm  x  17  cm  water  channel.  Various  arrangements 
of  roughness  elements,  shown  in  Fig.  1,  were  used  on  the  wall.  Velocity 
profiles  were  measured,  and  the  skin  friction  or  equivalent  sand  grain 
roughness  was  derived  from  the  logarithmic  portion  of  the  profile.  The 
simple  shapes  and  regular  spacing  of  the  roughness  elements  can  be  simu¬ 
lated  quite  well  by  our  roughness  model,  with  the  exception  of  the  short 
angles,  which  were  not  investigated. 

Figure  2  compares  the  skin  friction  computed  with  our  model  against 

Schlichting' s  data  for  the  spherical  segments.  As  indicated  in  Fig.  1, 

the  segments  are  nearly  hemispherical  (height  =  0.26  cm,  radius  =  0.40  cm). 

2 

Further,  the  dashed  lines  were  not  actually  presented  by  Schlichting.  He 
determined  the  equivalent  sand  grain  roughness,  based  on  the  observed  in¬ 
crease  in  Cf  over  the  smooth  wall  value  at  selected  downstream  stations. 

We  took  the  liberty  of  using  the  kg  values  to  obtain  the  augmentation  of 

as  a  function  of  distance.  In  so  doing,  the  appropriate  upstream  initial 
conditions  are  ignored.  The  computations  were  started  with  a  fully  turbulent 
smooth  wall  boundary  layer  at  x/Z  =  0.03.  Given  the  arbitrary  treatment  of 
initial  conditions,  the  mild  disagreement  at  upstream  locations  is  not  sig¬ 
nificant.  Otherwise,  the  model  reproduces  fairly  well  the  observed  increase 
in  skin  friction  with  increasing  roughness  density. 

A  detailed  comparison  of  the  mean  velocity  profiles  computed  from 
the  present  model  with  Schlichting' s  data  for  spherical  segments  is  shown 
in  Fig.  3,  in  semi-logarithmic  coordinates  (U^  =  A  /p^) .  Agreement  is 
again  seen  to  be  quite  good,  with  a  10-15%  error  for  the  most  dense  packing. 
The  data  and  curves  of  Fig.  3  agree  well  with  the  velocity  profile  of  classic 
rough  wall  pipe  flow  when  replotted  against  y/kg. 

Figure  4  shows  the  computed  skin  friction  coefficients  and  Schlich- 

ting’s  measurements  for  three  different  element  shapes  as  a  function  of  spac- 

7 

ing.  The  calculations  were  evaluated  at  Rex  s  10  ,  but  the  increase  in  skin 
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Fig.  3  Computed  mean  velocity  profiles  compared  to  Schlichting 
measurements  for  spherical  segments  at  three  relative 
spacings . 
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friction  is  essentially  independent  of  distance.  The  decrease  in  C  for 
spheres  as  2./D  *  1  is  easy  to  explain  qualitatively.  At  large  1/D  the  flow 
"sees”  the  entire  sphere.  With  substantial  element  spacing,  the  drag  in¬ 
creases  with  decreasing  Jl/D  because  there  are  more  elements  per  unit  area. 

But  as  £/D  -*■  1  the  flow  becomes  negligible  below  the  centers  of  the  spheres 
and  the  drag  is  due  only  to  the  upper  half  of  the  elements.  The  cones  and 
segments  were  not  investigated  at  such  close  spacings.  The  calculations 
for  the  cones  and  segments  are  almost  identical ,  although  the  cones  are 
slightly  taller,  0.375  cm  vs  0.26  cm.  However,  the  data  for  the  cones  fall 
above  that  for  the  segments.  This  implies  that  the  effective  drag  coeffi¬ 
cient  for  the  flow  about  a  conical  roughness  element  is  somewhat  larger  than 
that  for  flow  about  spherical  elements.  The  same  value  was  used  in  all  of 
our  calculations.  If  one  were  to  allow  such  a  higher  value  of  for  cones, 
the  influence  of  roughness  element  shape  would  be  well  understood. 

The  second  set  of  interesting  tests  was  carried  out  by  Acurex  Corpor¬ 
ation  in  AEDC  Tunnel  F,7  using  45°  conical  models  with  a  variety  of  surface 
roughnesses.  This  facility  was  an  arc-driven  hot  shot  tunnel,  in  which  the 
test  section  pressure  decreased  during  the  run  (total  time  =  200  msec) .  The 
most  useful  tests  were  performed  on  sharp  45°  cones  at  =  7  at  a  free  stream 
Reynolds  number  of  45  x  106/ft.  The  resulting  boundary  layer  edge  Mach  number 
was  1.7.  was  the  test  gas.  The  first  0.75"  of  wetted  length  was  roughened 
to  4-5  mils  to  ensure  rapid  transition. 

Seven  surface  finishes  were  used  on  the  remainder  of  the  cone-' 
essentially  smooth,  grit  blasted  to  almost  2  mils,  2  mil  bonded  grit,  and 
four  chemically- etched  roughness  patterns  (wide  and  close  spacing  at'  nominal 
heights  of  4  and  10  mils) .  The  etching  process  resulted  in  roughness  elements 
that  are  best  approximated  as  truncated  cones  (top  radius  =  1/4  base  radius) 
whereas  the  grit  roughness  elements  are  simulated  with  hemispheres.  The 

„  w 

roughness  characteristics  of  these  chemically  milled  surfaces  varied  by  as 
much  as  ±  30%  in  mean  roughness  height  over  the  surface  of  the  cones.  Table 
I  lists  the  average  element  height  and  spacing  for  the  various  surfaces.  Note 
that  the  "4  mil"  etched  roughness  is  actually  considerably  less  rough.  Also, 
the  "4  mil"  roughnesses  have  larger  relative  roughness  spacing  than  the  "10 
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TABLE  X. 


Roughness  Characteristics  for  Acurex  AEDC  Tests  at  Mach  7 


Designation 

Mean 
k  (mils) 

Mean  Spacing 

2.  (mils) 

.  V.  (-H 

\  m  sec/ 

Smooth 

— 

- ■ 

2038 

Grit  Blasted 

1.63 

7.70 

2235 

Bonded  Grit 

2.00 

4.00 

2292 

4  mil  Wide 

3.00 

23.0 

2157 

4  mil  Close 

2.50 

13.0 

2357 

10  mil  Wide 

10.32 

56.0 

1998  ' 

10  mil  Close 

9.50 

31.0 

2315 
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mil"  surfaces.  In  fact,  the  "4  mil  close"  and  "10  mil  wide"  roughnesses 
have  essentially  the  same  £/k. 

The  primary  measurements  for  this  test  series  are  the  heat  trans¬ 
fer  rates,  determined  by  thin  wall  calorimetry  with  thermocouples  on  the 
back  wall.  Some  co-axial  heat  transfer  gauges,  as  well  as  skin  friction 
gauges,  were  also  used,  although  there  may  be  uncertainties  regarding  how 
faithfully  the  roughness  is  reproduced  on  the  surface  of  the  gauges.  The 
skin  friction  measurements  are  limited  and  show  such  scatter  that  few  con¬ 
clusions  are  possible  (see  Ref.  8  for  more  discussion). 

Figure  5  compares  the  Stanton  numbers,  q/p  u  C  (T  -  T  ),  calculated 

^  e  e  p  r  w 

by  the  present  model  with  the  Acurex  data.  The  agreement  is  seen  to  be 
good,  with  the  theory  generally  well  within  the  scatter  of  the  data.  Several 
trends  are  evident  from  either  data  or  calculations.  The  bonded  grit  and 
grit  blasted  surfaces  cause  a  similar  heating  augmentation,  although  the  bonded 
grit  is  slightly  taller  and  considerably  more  densely  packed;  a  blockage  effect 
must  be  counter-acting  the  more  obvious  effect  of  element  spacing.  However, 
for  the  chemically  etched  surfaces,  spacing  appears  to  be  more  important  than 
height.  The  10-mil  wide  spacing  yields  a  greater  heating  rate  than  the  4-mil 
(actually  3-mil)  wide  case  only  at  larger  distances,  and  the  4-mil  and'  10-mil 
close  spacing  results  are  also  quite  similar.  As  already  noted,  the  relative 
roughness  spacing  is  not  constant  between  the  4-  and  10-mil  heights,  and  a 
more  detailed  discussion  of  the  dependence  of  height,  shape  and  spacing  is 
presented  below. 

Another  interesting  experiment  on  roughness  effects  has  been  con- 

9 

ducted  in  the  hypersonic  shock  tunnel  at  Calspan.  Only  one  roughness  was 
studied  —  a  bonded  grit  similar  to  that  employed  by  Acurex  in  the  series 
discussed  above.  The  mean  roughness  height  was  3.8  mils,  with  a  spacing 
of  10-15  mils  (we  specified  12.5  mils  for  the  spacing).  As  with  the  Acurex 

*The  data  points  have  actually  been  derived  by  dividing  the  reported  heat 
transfer  coefficients  by  our  computed  values  of  Peue  and  the  specific  heat 
of  N2-  Table  I  includes  the  values  used  for  peue  for  the  various  cases. 

Note  that  peue  varies  by  as  much  as  18%  from  one  case  to  another,  and  ex¬ 
amination  of  the  heat  transfer  rate  (q)  rather  than  the  Stanton  number 
could  possibly  lead  to  inaccurate  conclusions  on  the  effect  of  roughness. 
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series,  the  models  were  45°  cones.  The  free  stream  Mach  number  was  11-13, 
although  this  larger  value  has  little  effect  on  the  edge  Mach  number,  due 
to  the  large  cone  angle.  T^/T^  considerably  less  in  the  Calspan  condi¬ 
tions.  Heat  transfer  was  measured  with  thin  film  gauges  covered  with  the 
surface  roughness.  Figure  6  shows  our  comparison  with  Holden's  results, 
for  smooth  and  roughened  surfaces,  on  a  sharp  45°  cone  at  the  highest  Rey¬ 
nolds  number  tested.  The  smooth  wall  boundary  layer  was  naturally  turbu¬ 
lent  near  the  nose.  It  is  evident  that  the  model  yields  heating  rates 
about  15%  higher  than  measured,  even  for  the  smooth  wall.  The  measured 
smooth  wall  Stanton  number  is  consistent  with  values  predicted  from  ac¬ 
cepted  engineering  methods  such  as  that  of  van  Driest10  (St  s  2.56  x  10  3) . 
The  reason  for  this  discrepancy  is  not  clear,  although  it  may  be  associated 
with  numerical  inaccuracies  that  arise  with  large  density  differences  across 
the  boundary  layer.  In  any  case,  the  present  model  predicts  a  degree  of 
roughness  heating  augmentation  comparable  to  that  observed. 

The  first  high  Mach  number  tests  with  distributed  surface  roughness 
11 

were  conducted  by  Keel  on  5°  sharp  cones  in  Tunnel  No.  2  at  the  Naval 

Surface  Weapons  Center,  at  M  *  5  or  M  =4.77.  Sand  grains  were  uniformly 

applied  to  the  model  with  epoxy,  yielding  roughness  heights  of  either  ?3 

mils  or  43  mils.  Element  spacing  measurements  were  not  reported,  and  we 

assumed  a  value  (&/D  =  2.5)  which  is  typical  of  the  bonded  grit  surfaces 

7  9 

constructed  by  Acurex  and  Calspan.  Skin  friction  and  heat  transfer  were 
measured  with  floating  element  balances  and  slug  calorimeter  gauges,  respec¬ 
tively. 

One  noteworthy  aspect  of  Keel's  experiment  is  that  the  measurements 
were  obtained  at  a  fixed  station  on  the  cone,  x  =  2  ft.  The  Reynolds  num¬ 
ber  was  varied  by  decreasing  the  tunnel  pressure,  and  the  results  plotted 

as  C.  or  St  vs  Rea  (see  Fig.  7) .  The  manner  by  which  the  data  were  collected 
f  o 

must  be  recognized  to  properly  interpret  these  results.  With  a  smooth  wall, 

C.  is  a  function  only  of  Ren  (assuming  fixed  values  of  M  and  T  /T  )  and  one 
f  0  ewe 

may  equally  well  traverse  the  Cf  vs  Reg  curve  by  varying  pressure  (i.e.,  den¬ 
sity)  or  distance  (i.e.,  6).  However,  a  second  independent  parameter  enters 
with  rough  walls.  For  simplicity,  let  us  use 


P  u  k/y 

W  T  W 
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rough  (4  mil  bonded  grit) . 


as  the  fundamental  parameter  to  describe  the  roughness  augmentation  effect. 

Furthermore,  u  /u  is  a  very  weak  function  of  Re.,  and  may  be  considered 
T  e  o 

constant  for  the  purposes  of  this  argument.  Thus,  if  we  consider  increas- 

0  8 

ing  distance  along  a  flat  plate  or  cone,  Reg  increases  because  0  -  x  '  , 
while  k+  remains  constant.  Conversely,  when  Keel  raises  the  pressure  at  a 
fixed  station,  both  Reg  and  k+  increase  together,  linearly  with  the  density. 
This  offers  an  explanation  for  the  almost  complete  absence  of  slope  for  the 
data  in  Fig.  7;  the  increasing  roughness  augmentation  (with  increasing  k+) 
tends  to  cancel  the  natural  tendency  of  or  St  to  decrease  with  increas¬ 
ing  Reg.  Agreement  with  the  data  is  very  good  -  k+  is  apparently  suffi¬ 
ciently  large  that  the  roughness  augmentation  is  nearly  saturated  and  there 
is  little  difference  between  23  and  43  mils.  It  should  be  noted  that  the 
computed  Mach  number  at  the  tops  of  the  elements  is  barely  supersonic  in 
this  case. 

For  hypersonic  edge  conditions.  Fig.  8  shows  a  comparison  with  the 

12 

heat  transfer  data  obtained  by  Hill  at  NSWC  on  a  7°  sharp  cone  at  M  =8 

6  ^ 
and  Rem  =  7  x  10  /ft.  The  roughness  heights  are  based  on  the  nominal  sizes 

of  the  grits  that  were  applied  to  the  surface.  We  simulated  the  roughness 

with  hemispheres  of  height  equal  to  the  nominal  values,  with  a  spacing  of 

2.9k,  as  indicated  by  surface  profilometer  measurements.  It  is  evident 

from  Fig.  8  that  the  computer  model  yields  approximately  the  observed  heat 

transfer  augmentation  as  a  function  of  roughness  height,  although  the  data 

show  less  dependence  on  distance.  The  reasons  for  this  discrepancy  are  not 

clear. 

An  even  more  intriguing  set  of  hypersonic  data  have  been  obtained  by 
Holden1^  in  the  Calspan  shock. tunnel .  In  this  case  the  edge  Mach  number  is 
9.4,  the  cone  half-angle  is  6°,  and  the  free-stream  Reynolds  number  for  the 
cases  of  interest  is  1.1  x  10^/ft.  Again,  we  simulated  the  roughness  with 
hemispheres  with  height  equal  to  the  nominal  grit  sizes.  In  this  case, 
however,  the  elements  are  rather  tightly  packed,  £/k  »  2.25.  Figure  9  shows 
that  the  theory  yields  a  trend  of  increasing  heat  transfer  with  increasing 
roughness,  while  the  data  show  significant  decreases.  There  are  several 
possible  explanations  for  the  observed  reductions  in  heat  transfer.  One 


St  =  q/pe  Ug  (H 


DISTANCE  ALONG  CONE  (INCHES) 


Fig.  8  Comparison  of  present  theory  with  Hill's  data  for 
smooth  and  rough  surfaces  at  an  edge  Mach  number 
of  8. 
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Mach  number  of  9.4. 


is  associated  with  the  tight  packing  of  the  elements.  However,  our  calcu¬ 
lations  show  that  there  might  be  a  modest  reduction  in  the  heating  augmen¬ 
tation  as  £/D  -*■  1  (cf.  Fig.  4),  but  not  below  smooth  wall  values.  Further¬ 
more,  Holden's  data  on  45°  cones  at  M  =  1.8,  with  a  similar  4  mil  bonded 

e 

grit  roughness,  were  presented  above  in  Fig.  6  and  showed  no  reduction 
below  smooth  wall  values.  A  second  possibility  has  to  do  with  upstream 
beh'avior;  the  rough  wall  boundary  layer  may  simply  be  much  thicker  at  a 
given  station.  Again,  the  present  calculations  do  not  indicate  any  such 
reversal  effect.  A  third  possibility  has  to  do  with  the  high  Mach  number, 
which  could  result  in  locally  supersonic  flow,  shock  waves  around  the  ele¬ 
ments,  and  perhaps  an  alteration  of  the  drag  or  heating  laws.  In  fact, 
Holden's13  shadowgraph  observations  indicate  a  significant  wave  structure 
emanating  from  the  rough  surface.  But,  our  computed  local  Mach  numbers 
are  only  barely  supersonic  (1.14)  at  the  tops  of  the  elements,  no  more 
supersonic  than  in  the  Keel  case  discussed  above.  A  related  Mach  number 
effect,  discussed  in  more  detail  below,  is  that  roughness  causes  increased 
boundary  layer  temperatures  at  high  edge  Mach  numbers;  perhaps  more  careful 
calculations  would  show  a  reduction  below  smooth  wall  values.  In  any  case, 
the  matter  is  not  resolved,  particularly  in  light  of  the  similar  conditions 
of  the  experiments  of  Hill  and  Holden.  Clearly,  more  measurements  and 
calculations  are  needed. 
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IV.  ROUGHNESS  SCALING  BEHAVIOR 


A  major  advantage  of  the  present  type  of  numerical  study  is  that 
the  solutions  can  be  scrutinized  to  determine  the  nature  of  roughness  ef¬ 
fects.  Of  course,  the  data  comparisons  shown  in  the  previous  section  left 
several  issues  unresolved,  but,  nevertheless,  a  detailed  examination  of 
the'numerical  results  is  justified.  The  results  are  quite  interesting. 

The  most  striking  aspect  of  our  solutions  is  that  the  mean  velo¬ 
city  is  computed  to  be  quite  uniform  over  much  of  the  range  y  <  k.  A 
typical  example  is  shown  in  Fig.  10.  Of  course,  very  near  the  wall  (i.e. , 
at  the  bottom  of  the  elements)  the  velocity  must  be  zero,  and  near  the 
tops  of  the  elements  the  velocity  tends  to  increase.  This  constant  velo¬ 
city,  which  was  unexpected,  is  evident  in  almost  all  of  the  cases  consider¬ 
ed,  the  sole  exceptions  being  cases  with  very  small  roughness  (say  k+  <  10) 
or  very  large  roughness  (k/9  >  1).  Within  the  range  of  constant  velocity, 
production  of  turbulent  energy  is  negligible;  the  turbulence  simply  diffuses 
toward  the  wall  and  dissipates.  In  the  mean  velocity  equation  (Eq.  3)  only 
the  final  two  terms,  representing  turbulent  shear  and  form  drag  on  the  ele¬ 
ments,  are  important.  The  total  drag  on  the  elements,  which  should  be 
close  to  the  actual  rough  wall  skin  friction,  is  easily  related  to  the 
element  shape  and  spacing 


Here  we  use  (  )  to  denote  conditions  in  the  region  of  constant  velocity. 

H 

The  integral  in  Eq.  6  is  the  frontal  area  of  the  elements  per  unit 

superficial  surface  area.  This  quantity  is  easily  related  to  the  shape  and 

density  of  roughness  elements.  However,  to  compute  the  skin  friction  the 

roughness  velocity  U  is  also  needed  (determination  of  p  ,  as  well  as  the 

R  K 

relation  between  heat  transfer  and  skin  friction,  will  be  addressed  below) . 
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Calculated  mean  velocity  profile  for  a  typical  rough  wall  case. 


We  have  been  unable  to  predict  U  by  any  simple  theoretical  arguments. 

Accordingly,  the  computer  model  was  systematically  run  over  a  range  of 

parameters,  including  element  shape,  spacing,  height,  T  /T  ,  and  M  .  To 

W  6  6 

be  specific,  U  was  defined  to  be  the  mean  velocity  at  y  =  k/2. 

K 

Figure  11  shows  the  computed  roughness  velocity  for  hemispheres 
in  incompressible  conditions.  The  roughness  spacing  varies  along  the  ab¬ 
scissa;  U  is  normalized  by  the  smooth  wall  velocity  at  y  =  k/2  (easily 
computed  from  the  classic  law-of-the-wall) .  Only  a  slight  dependence  on 
Reg  is  apparent.  It  should  be  emphasized  that  only  roughness  spacing  and 
downstream  distance  (Reg)  were  varied  to  obtain  these  curves.  However, 
the  really  significant  finding  is  that  the  same  curves  apply  to  variations 
in  element  shape,  height,  T^/Tg,  and  M^.  In  addition  to  hemispheres,  we 
also  considered  30°  conical  elements  (base  diam  =  k) ,  45°  cones  (base  diam 
*  2k),  truncated  30°  cones  (base  diam  =  k,  top  diam  =  k/2),  and  cylinders 
(base  diam  =  k) .  For  each  of  these  shapes  the  results  were  indistinguish¬ 
able  from  those  of  Fig.  11.  There  is  a  dependence  on  roughness  height  (k+) , 
as  illustrated  in  Fig.  12,  but  the  trend  is  slight  unless  k+  is  less  than 

about  20  (as  k+  -*•  0,  U  /U  ■*  1,  but  the  smooth  wall  shear  would  dominate 

R  S 

anyway).  Varying  the  wall  temperature  over  the  range  T /T  =  0.2  -  1.0  has 

W  6 

little  effect,  as  does  increasing  the  edge  Mach  number  up  to  at  least  8. 

The  substantial  effect  of  roughness  on  the  velocity  profile  is  not 
observed  with  the  enthalpy  profile.  This  undoubtedly  results  from  the 
fact  that  there  is  no  heat  transfer  mechanism  analogous  to  form  drag.  In 
Fig.  13  we  compare  normalized  total  enthalpy  values  at  y  =  k/2  with  the  rough¬ 
ness  velocity.  At  =0,  the  rough  and  smooth  wall  enthalpies  at  y  =  k/2 
are  identical  within  the  accuracy  of  our  computations,  even  for  large  tem¬ 
perature  differences  across  the  boundary  layer  (T  =  T  /5).  At  supersonic 

W  0 

velocities,  the  total  enthalpy  is  obviously  reduced  by  smaller  margins  than 

-  ^  ^  1/2 

is  the  velocity.  There  appears  to  be  a  hypersonic  limit,  with  H  /H  ~  (u  /U  )  , 

R  5  R  S 

although  more  computations  would  be  required  to  confirm  this. 
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PROJECTED  ELEMENT  AREA  /  UNIT  AREA 


Ratio  of  computed  roughness  velocity  (y  =  k/2)  to  smooth  wall  velocity  at  the  same 
height,  as  a  function  of  projected  roughness  element  area  and  k+. 


Ratio  of  computed  total  enthalpy  (y  =  k/2)  to  smooth  wall  enthalpy  at  the  same 
height,  as  a  function  of  projected  element  area  and  edge  Mach  number. 


An  important  consequence  of  Fig.  13  is  that  roughness  will  tend  to 
increase  the  static  enthalpy  for  y<k  at  supersonic  velocities.  The  resulting 
decrease  in  density  will  reduce  the  form  drag  on  the  elements,  and  the  in¬ 
creased  temperature  will  reduce  the  underlying  smooth  wall  skin  friction. 

Thus  our  calculations  suggest  a  specific  Mach  number  dependence  for  rough¬ 
ness  .friction  or  heating  augmentation.  The  compensating  effect  of  increased 
temperatures  reduces  roughness  augmentation  at  high  Mach  numbers;  if  this  is 
somewhat  greater  than  we  have  estimated,  it  could  offer  an  explanation  for 
the  decreases  in  heat  transfer  observed  by  Holden. 

Finally,  we  might  note  two  additional  trends  that  are  well  established 
by  our  numerical  results.  First,  the  underlying  smooth  wall  skin  friction 
is,  to  a  good  approximation,  unchanged  by  roughness  (fixed  value  of  Reg) » 
except  for  the  high  Mach  number  effect  just  mentioned.  Secondly,  the  increase 
in  heat  transfer  is  nicely  described  by  the  square  root  of  the  skin  friction 
augmentation.  The  reason  was  given  in  our  previous  study^:  roughness  in¬ 
creases  fluctuating  velocities  but  not  fluctuating  temperatures,  and  the 
turbulent  stress  goes  as  -pu’v'  while  the  turbulent  heat  flux  goes  as  -pu’T’  . 


a 
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V.  ROUGHNESS  SCALING  FORMULAE 

It  is  useful  to  summarize  the  relations  needed  to  predict  roughness 
effects  on  skin  friction  and  heat  transfer,  particularly  since  the  present 
results  are  in  a  very  different  form  from  existing  methods.  First,  the 
skin  friccion  and  heat  transfer  augmentations  are  given  by 


„  „  PR  U: 

Cf=Cf  +— 

S  P_ 


r  r  f  (k/2) 

u  i  a 


(7) 


(8) 


The  smooth  wall  values  may  be  obtained  from  a  variety  of  well-known  sources, 
such  as  Ref.  10.  The  integral  in  Eq.  (7)  requires  the  roughness  element 
shape  and  spacing,  and  f(k/2)  follows  from  Eq.  (4).  We  recommend  CD  =  0.6. 
The  ratio  UR/Ug  may  be  obtained  from  Figs.  11  or  12,  and  the  smooth  wall 
velocity  at  y  =  k/2  from  the  standard  law-of- the -wall : 


(9) 


The  roughness  density  appearing  in  Eq.  (7)  requires  slightly  more 
algebra.  Assuming  a  perfect  gas, 
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It  is  adequate  to  use  a  Reynolds  analogy  for  the  smooth  wall  total  enthalpy, 

(H-h  )/(H~-h  )  =  U/U  ,  to  estimate  the  smooth  wall  static  enthalpy: 
w  u  w  e 
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_§  = + 
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And,  the  total  enthalpy  ratio  in  Eq.  (10)  is  directly  related  to  the  nor¬ 
malized  total  enthalpy  difference 
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Taking  H  /H  from  Fig.  13,  Eqs.  (10)- (13)  are  sufficient  to  provide  P  /P  . 

R  S  R  6 

Finally,  at  high  Mach  numbers,  the  temperature  increase  due  to 
roughness  can  reduce  the  smooth  wall  friction  coefficient  or  Stanton  num¬ 
ber.  To  estimate  this  effect,  we  suggest  modifying  a  simple  reference 
enthalpy  method.  The  smooth  wall  coefficients  are  inversely  proportional 
to  the  reference  enthalpy,  which  in  turn  depends  on  the  peak  boundary  layer 
enthalpy 


ref 


s  1  +  0 


(14) 


For  a  smooth  wall,  h  may  be  estimated  by  taking  the  maximum  of  Eq. 
(12)  with  respect  to  U  :  - 
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However,  if  the  rough- wall  enthalpy  resulting  from  Eq.  (10)  is  greater 
tan  the  value  obtained  from  Eq.  (15),  then  the  smooth  wall  friction  co¬ 
efficient  and  Stanton  number  should  be  reduced  by  the  factor 


i  ♦  -  i 

— (H  ' 


(16) 
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VI.  DISCUSSION 


The  present  rough  wall  Reynolds  stress  model  yields  results  that 
compare  acceptably  with  a  variety  of  available  measurements.  Of  the  sev¬ 
eral  discrepancies  noted,  the  most  significant  is  that  with  Holden's  high 
Mach  number  tests.  The  matter  may  not  be  resolved  until  additional  high 
Mach  number  data  are  obtained.  The  observed  decrease  in  heating  rates 
could  be  due  to  the  extreme  close  packing  of  the  roughness  elements,  al¬ 
though  our  model  fails  to  confirm  this.  The  model  does  suggest  that  the 
increased  temperatures  generated  by  roughness  at  high  edge  Mach  numbers 
tend  to  reduce  skin  friction  and  heat  transfer.  However,  these  findings 
are  only  suggestive,  and  one  cannot  at  this  time  rule  out  an  experimental 
artifact  or  a  significant  change  in  the  nature  of  flow  about  .roughness 
elements  in  supersonic  conditions. 

The  roughness  scaling  behavior  derived  from  our  systematic  vari¬ 
ation  of  the  several  relevant  parameters  suggests  scaling  laws  that  are 
a  major  departure  from  previous  approaches.  The  influence  of  roughness 
element  shape  and  spacing  is  almost  totally  defined  by  the  projected  front¬ 
al  area  per  unit  surface  area.  The  velocity  that  gives  rise  to  form  drag 
is  largely  constant  over  the  y  <_  k.  When  normalized  by  the  smooth  wall 
velocity  at  y  *  k/2,  this  roughness  velocity  is  primarily  a  function  of 
this  projected  element  area/unit  area. 

A  simple  method  is  thus  provided  for  predicting  and  analyzing  rough¬ 
ness  character  effects.  However,  it  must  be  admitted  that  the  relations 
given  here  need  additional  evaluation  and  refinement  before  they  can  be 

used  with  great  confidence.  For  example,  the  minor  dependence  of  U  /U 

+  R  s 
on  k  and  Reg  should  be  better  defined  by  additional  calculations,  and 

much  more  hypersonic  data  ought  to  be  studied  to  substantiate  the  Mach 

number  dependence. 
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APPENDIX 


REYNOLDS  STRESS  MODEL  EQUATIONS 


The  governing  partial  differential  equations  for  the  various  tur¬ 
bulent  quantities  ere  listed  here.  For  more  details,  the  reader  is  re- 

ferred  to  our  previous  report. ^  In  practice ,  it  is  convenient  to  replace 

2  2  2  2  2  2  2 
u*  ,  v'  ,  w1  by  the  kinetic  energy  q  =  (u1 _ +  v'  +  w'  ) /2  and  two 

2  2  2  2 
measures  of  the  degree  of  anisotropy  S^  =  u'  -2/3  q  ,  S 22  =  v‘  -  2/3  q 

In  all  of  the  following  equations,  the  diffusion  terms  (second  partial 

derivatives  with  respect  to  y)  should  be  understood  to  be  multiplied  and 

divided  by  B(y)  as  in  Eq.  3. 

For  steady  flow,  the  governing  equations  include  continuity: 

it  (pUi>  •  0  u-i, 


the  mean  momentum  equation: 
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1.  2  +  12.  5  TT/Re^ 
where  C£  =  j  +  l2.  5  TT/ReA  "  ' 

(o.  288  +  6.  6  TT/ReA  +  35  TT2/1^2) 
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5  (0.4  +  5  TT/Re^r 

and  Re^  is  the  turbulent  Reynolds  number  qA/v,  with  A  being  related  to  the 
dissipation  rate  by 


37 


0.4 


+  5  rrv 


0.  4  (1  +  12.  5  TT/Re.  ) 

A  A 


S 


the  corresponding  equations  for  enthalpy- related  quantaties  are: 
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The  terms  R^,  R^_,  R^,  and  R^  contain  the  effect  of  roughness  on 
the  boundary  layer.  For  the  mean  velocity  and  enthalpy  equations,  R^ 
and  Rh  were  presented  in  Eqs.  (1)  and  (2)  in  Section  II,  the  other  two 
terms  are  source  terms  for  kinetic  energy  and  dissipation,  describing  the 
fluctuations  introduced  in  the  wakes  of  elements.  For  the  fully  turbu¬ 
lent  boundary  layers  considered  in  this  study,  these  terms  are  generally 
small  compared  to  the  natural  turbulence  production  terms.  As  described 
in  Ref.  5,  the  terms  used  are: 


R  =  0.04pu3D/il2 
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